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This  is  the  Final  Report  to  ARO  Grant 
"Theory  of  Heterojunction  Discontinuities" 

November  1976  - 31  December  1977.  The  Principal  Investigator 
was  Dr.  Herbert  Kroemer,  Professor  of  Electrical  Engineering  at 
the  University  of  California  at  Santa  Barbara  (UCSB) . The 
grant  was  a follow-up  to  an  earlier  grant  (#  DAHCO4-74-G-0114) 
to  the  University  of  Colorado,  with  the  same  Principal  Inves- 
tigator. A Final  Report  on  the  predecessor  grant  was  issued  in 
December  1976;  the  present  report  covers  the  follow-up  grant 
only. 
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2 . ) The_Problem 


The  problem  studied  was  the  theoretical  understanding  of 
the  energy  band  lineup  at  semiconductor  hetero junctions  and  the 
development  of  theoretical  means  to  predict  that  lineup  for  new 
and  as  yet  untried  hetero junction  pairs. 

Hetero junctions  are  junctions  between  different  semicon- 
ductors, such  as  GaAs  and  AlAs.  They  are  increasingly  being 
employed  in  advanced  semiconductor  devices,  and  the  performance 
of  the  devices  depends  on  the  relative  lineup  of  the  energy 
bands  at  the  junctions.  Prior  to  the  work  under  the  present 
grant  and  its  predecessor,  the  only  rule  for  the  prediction  of 
the  lineup  was  the  electron  affinity  rule.  It  is  theoretically 
ill-founded,  and  in  practice  hard  to  apply,  because  it  requires 
empirical  surface  property  data  on  high-quality  crystals.  Even 
where  such  data  are  available,  the  predictions  made  from  them 
often  do  not  agree  with  the  results  obtained  on  actual  hetero- 
junctions — once  the  latter  have  been  prepared,  often  at  great 
technological  expense. 

What  was  desired  was  a better  theoretical  understanding  of 
why  the  hetero junctions  line  up  the  way  they  do,  on  a level 
that  would  make  it  possible  to  predict  the  lineup,  with  a 
minimum  of  technological  investment,  and  thus  to  aid  in  making 
decisions  as  to  whether  or  not  the  odds  are  favorable  to  make 
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it  vrorthwhile  to  embark  on  the  possibly  very  expensive  technology 
to  build  a particular  new  device.  We  succeeded  in  developing 
theoretical  techniques  that  have  such  a predictive  capability. 

Under  the  predecessor  grant  this  work  had  reached  the 
following  stage. 

A self-consistent  pseudopotential  method  was  developed 
that  attempted  to  account  for  the  known  ionic  core  potential  as 
realistically  as  possible  before  applying  a repulsive  pseudo- 
core to  the  overall  potential.  By  matching  the  calculated 
energy  band  structure  of  each  semiconductor  to  the  observed 
band  structure,  it  was  possible  to  determine  the  electrostatic 
potential  outside  the  cores  that  goes  hand-in-hand  with  this 
band  structure.  This  permits  expressing  the  band  structure 
relative  to  a suitably  chosen  reference  point  on  the  electro- 
static potential,  a key  step  in  determining  hetero junction 
lineups.  Such  a reference  potential  should  be  chosen  as  far 
away  from  the  ionic  cores  as  possible.  Two  kinds  of  points 
suggested  themselves:  the  halfway  point  along  the  interatomic 
bonds,  and  the  points  in  the  centers  of  the  large  interstitial 
sites.  There  are  two  kinds  of  interstitial  sites  in  the  compounds 
of  interest,  suggesting  the  average  of  the  two  interstitial 
potentials  as  a reference  potential.  When  we  made  the  ad  hoc 
assumption  that  the  latter  reference  potential  be  continuous 
across  the  heterojunction,  we  obtained  lineups  that  agreed 
surprisingly  well  with  observed  lineups,  much  better  than  for 
other  choices  of  the  reference  potential.  However,  the  choice 
seemed  to  lack  a good  physical  foundation.  The  best  justifi- 
cation we  could  give  it  was  a muffin-tin  potential  analogy  that 
was  not  very  satisfactory.  Our  first  objective  for  the  follow- 
up grant  was  to  come  up  with  a more  rational  electrostatic 
potential  matching  scheme. 

3 . ) 

We  discovered  early  in  the  present  grant  that  the  mean 
interstitial  potential  is  the  potential  at  infinity  if  one 
treats  the  semiconductor  as  an  assembly  of  spherical  ions 
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bounded  only  by  non-polar  surfaces,  and  that  it  is  therefore 
indeed  the  natural  reference  potential.  Effective  ionic  charges 
for  various  semiconductors  were  determined  and  were  found  to 
correlate  strongly  with  the  difference  in  Phillips  electronega- 
tivities between  the  constituent  atoms  of  the  semiconductor 
[1] . If  one  allocates  one-quarter  of  each  ionic  charge  to  a 
charge  transfer  along  each  bond  one  can  associate  a charge 
transfer  with  each  bond.  If  one  applies  this  charge  transfer 
to  the  bonds  crossing  the  interface  one  obtains  an  estimate  for 
an  interface  dipole  by  which  the  reference  potentials  on  the 
two  sides  are  shifted  relative  to  each  other  [2]. 

Towards  the  end  of  the  grant  we  modified  the  spherical  ion 
model  by  assuming  that  there  are  actually  charges  located  along 
the  bonds.  This  leads  to  small  corrections  in  the  lineups,  and 
to  a small  orientation-dependence  of  the  hetero junction  dis- 
continuity. The  correction  is  zero  for  (110)  planes,  has 
opposite  values  for  (111) A and  (lll)B  planes.  The  values  for 
the  two  (100)  planes  are  the  same  in  this  model  as  for  the 
corresponding  (111)  planes.  Numerical  values  can  not  be  esti- 
mated at  this  point  because  of  the  absence  of  accurate  values 
for  the  magnitude  and  location  of  the  bond  charges.  This  last 
work  has  not  been  published  yet;  a preprint  is  attached  as 
Appendix. 

During  the  last  three  months  of  the  grant  we  briefly 
considered  interface  reconstruction  effects  and  effects  of 
grading,  as  in  LPE  junctions.  No  corrections  of  any  signifi- 
cance were  found. 
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APPENDIX 

A THEORETICAL  STUDY  OF  THE  ORIENTATION  DEPENDENCE 
OF  HETEROJUNCTION  ENERGY  BAND  LINEUPS* 

By 

VI.  R.  Frensley 
Texas  Instruments  Inc. 

Dallas,  Texas  75222 

l,i 

II  ABSTRACT 

Li 

{* 

\ The  orientation  dependence  of  the  energy  band  discontinuities  at  an  abrupt 

: semiconductor  heterojunct ion  is  studied  by  modeling  the  microscopic  charge 

l]  distribution  near  the  heterojunction.  The  model  includes  both  ionic  and  bond 

tj  charges,  and  assumes  that  each  chemical  bond  contributes  independently  to  the 

total  charge  distribution.  The  energy  band  lineups  are  derived  by  calculating 
the  interface  dipole.  The  results  indicate  that  there  should  be  no  difference 
in  the  band  lineups  for  heterojunctions  on  (100)  and  (111)  faces  of  similar 
polarity.  The  lineup  should  be  completely  independent  of  orientation  for  heterojunc- 
tions  which  have  an  element  common  to  both  semiconductors,  such  as  GaAs-AlAs.  Any 
difference  in  band  lineup  between  heterojunctions  on  nonpolar  (110)  planes 
and  those  on  the  (100)  or  (111)  planes  is  related  to  deviations  of  the  bond 
parameters  from  chemically  systematic  behavior,  and  thus  is  expected  to  be  small. 

PACS  numbers:  73-^0  Lq,  73«30  + y 
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I . I ntr  diict ion 

An  important  property  of  an  abrupt  semiconductor  he teroj unction  is 
the  change  in  energy  band  gap  across  the  interface,  which  leads  to  dis- 
continuities in  the  band-edge  energies.  The  band  gap  change  is  distributed 
between  the  conduction  and  valence  band  discontinuities.  The  resulting  band 
lineup  is  a property  of  the  heterojunction,  and  it  is  an  important  consideration 
in  the  choice  of  heterojunction  materials  for  a given  device  application. 

We  have  previously  proposed  a method  for  theoretically  predicting  the 
band  lineup  at  an  abrupt  heterojunct ion. ^ We  divided  the  lineup  problem  into 
two  parts.  The  first  part  concerned  the  calculation  of  the  electronic  band 
structure  of  the  bulk  semiconductors,  and  the  relation  of  the  band  structure  j 

to  the  electrostatic  crystal  potential.  The  second  part  concerned  the  lineup 
of  the  electrostatic  potential  across  the  heterojunction,  which  was  calculated 
from  a model  of  the  microscopic  charge  distribution  near  the  interface.  In 
the  present  work,  we  will  investigate  a more  realistic  model  of  the  charge 
distribution,  and  emphasize  the  differences  in  band  lineup  which  might 
occur  between  heterojunctions  of  similar  composition,  but  which  are  fabricated 
on  different  crystallographic  planes.  Note  that,  within  the  above  division  of 
the  problem,  any  orientation  dependence  must  be  due  to  differences  in  the  charge 
distribution,  not  the  bulk  band  structure.  ' 

r II.  Independent  Bond  Model 

In  Reference  1 we  took  the  valance  charge  distribution  to  be  a linear 
superposition  of  spherical  charges  centered  on  the  atomic  sites.  Such  a 

model  neglects  an  important  feature  of  the  true  charge  distribution  in  ^ 
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tetrahedral  semiconductors,  which  is  the  accumulation  of  charge  between  the 
atoms.  Brill  took  this  accumulation  of  charge  into  account  in  his  X-ray  analysis  of 
diamond  by  introducing  spherical  bond  charges  midway  between  the  carbon 
atoms, ^ and  Phillips  emphasized  the  structural  role  of  the  bond  charge.^ 

Walter  and  Cohen  calculated  the  valence  charge  density  of  seven  semiconductors 
using  wavef unct ions  determined  by  the  empirical  pseudopotential  method.'* 

Their  results  indicate  that  the  charge  distribution  can  be  well  approximated 
by  a superposition  of  spherical  ionic  and  bond  charges.  As  the  ionicity 
of  the  crystal  increases,  the  center  of  the  bond  charge  moves  toward  the 
anion  and  the  magnitude  of  the  bond  charge  decreases.  Of  course,  as  the 
ionicity  increases  more  of  the  valence  charge  is  found  in  the  vicinity 
*of  the  anion. 

Recently,  calculational  techniques  have  been  developed  which  permit  a 
self-consistent  evaluation  of  the  valence  charge  distribution  in  the  vicinity 
of  the  surfaces^'®  and  interfaces.^'®  The  results  of  these  calculations 
show  some  striking  features.  First,  the  perturbation  in  the  charge  distribution 
due  to  an  interface  extends  over  only  a very  short  distance.  Near  a surface. 


the  charge  distribution  assumes  the  bulk  characteristics  by  the  third  atomic 
layer.®*®  At  a heterojunction  the  disturbance  is  even  more  short-ranged.  Baraff, 
Appelbaum,  and  Hamann^  have  performed  detailed  calculations  on  a system 
consisting  of  three  atomic  layers  of  Ge  on  a semi - i nf i ni te  GaAs  crystal.  They 
found  that  the  charge  distribution  in  the  vicinity  of  the  second  Ge  layer  was 


indistinguishable  from  that  of  bulk  Ge.  Also,  the  bond  charges  between  the  Ge 
and  Ga  atoms  had  a magnitude  intermediate  between  those  of  bulk  Ge  and  GaAs. 
Moreover,  they  were  able  to  determine  the  effective  ionic  charge  associated  with  the 
interface  bonds,  and  found  that  it  was  equal  to  one-half  the  corresponding  value 


for  bulk  GaAs.  All  of  this  lends  support  to  the  notion  that  one  can  regard  the 
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total  valence  charge  distribution  as  a superposition  of  bond  distributions, 
and  that  the  chemical  trends  observed  in  the  bulk  semiconductors  also  hold 
for  bonds  at  a heterojunction  interface. 

In  this  paper,  we  will  conjecture  that  such  behavior  is  a characteristic 
of  all  heterojunctions.  Specifically,  we  will  assume  the  following: 

(i)  the  microscopic  charge  distribution  is  a linear  superposition  of 
individual  charge  distributions,  each  of  which  is  associated  with  a single 
chemical  bond,  and  (ii)  those  charge  distributions  are  functions  only  of  the 
chemical  species  being  bonded.  Thus,  we  assume  that  each  bond  contributes 
independently  to  the  total  valence  charge  distribution.  We  should  emphasize 
that  this  is  different  from  the  approximations  made  in  the  linear-combination- 
of-atomi c-orbi ta 1 s (LCAO)  theory  of  heterojunction  lineups.^  We  do  not 
assume  that  different  bonds  are  not  coupled  by  Hamiltonian  matrix  elements, 
but  we  do  assume  that  the  complete  self-consistent  solution  to  the  interface 
charge  density  can  be  resolved  into  individual  bonds. 

The  specific  model  of  the  microscopic  charge  distribution  which 
we  wish  to  investigate  is  described  as  follows.  The  charge  distribution 
is  assumed  to  be  a linear  superposition  of  spherical  ionic  and  bond 
charges,  with  the  ionic  charges  centered  on  the  equilibrium  atomic  sites 
and  the  bond  charges  centered  on  some  point  along  the  line  joining  nearest 
neighbor  atomic  sites.  The  position  of  the  center  of  the  bond  charge, 
measured  from  the  anion,  is  denoted  x and  is  expressed  as  a fraction  of 
the  bond  length  /j  a/^,  where  a is  the  cubic  lattice  constant.  Thus,  the  microscopic 
charge  distribution  may  be  completely  described  by  the  total  charge  contained 
in  the  ionic  and  bond  charges,  the  width  of  their  respective  charge  distributions, 
and  the  bond  charge  position  x.  If,  however,  we  are  only  interested  in  evaluating 
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the  electrostatic  potential  at  a charge-free  point  (such  as  the  interstitial 
point),  we  may,  without  loss  of  generality,  take  the  width  of  the  charge 
distributions  to  be  zero  (that  is,  to  assume  point  charges).  Therefore, 
we  need  only  to  consider  the  magnitude  of  the  ionic  and  bond  charges  and 
the  position  of  the  bond  charge. 

Since  the  crystal  must  be  electrically  neutral  and  each  bond  contributes 
independently  to  the  charge  distribution,  the  charge  distribution  attributed 
to  each  bond  must  be  neutral  overall.  This  distribution  consists  of 
three  charges  (anion,  cation,  and  bond  charges),  so  the  magnitudes  of  those 
charges  must  depend  on  only  two  parameters.  The  calculations  will  be 
simplified  if  we  choose  as  parameters,  the  bond  charge,  denoted  q.  and  a 

D 

transferred  charge,  denoted  q^.  We  assume  that  the  bond  charge  is 
derived  from  both  the  anion  and  cation,  with  a charge  xq„  coming  from 

D 

the  cation  and  (l-x)q„  coming  from  the  anion.  The  charge  q_  is  transferred 
from  cation  to  anion.  Thus  the  charge  on  the  anion  due  to  a single  bond  is 


"anion  ’ "t  ' "b’ 


and  that  the  cation  is 


(1) 
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‘’cation  “ "‘’t  ' (2) 

as  shown  in  Figure  1.  The  total  charge  on  a given  ion  is  simply  the  sum 
of  the  charges  due  to  each  of  the  four  chemical  bonds  involving  that  ion. 

We  have  previously’  defined  the  dipole  potential  as  the  difference  j 

between  the  mean  interstitial  potentials  V,  on  either  side  of  the  heterojunction.  j 


Vp  - V. (1)  - V.  (2). 
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For  the  present  calculations,  it  will  be  I’ore  convenient  to  consider  the 
spatially  averaged  potential  V and  its  difference  across  the  heterojunction 

AV  = V(l)  - V(2) 

The  quantities  V.  and  V can  bt  ''e!ated,  within  the  bulk  semiconductor,  by  a 
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Hide  1 ung- type  calculation.  We  have  previously  described  such  a calculation 
for  a spherical  ion  model. To  take  the  bond  charges  into  account,  we  consider 
a lattice  of  charges  q-  located  at  the  bond  charge  sites  and  charges  ■2q- 

D D 

located  at  the  atomic  sites.  The  potential  at  an  interstitial  point  of  this 
lattice  can  be  expressed  by 


6.  (x)  q„ 
o 


(5) 


where  i = 1,2  denotes  the  two  inequivalent  intersticies.  The  values  of  6 
are  evaluated  with  the  boundary  condition  that  the  spatially  averaged  potential 
equals  zero.  Therefore,  we  can  relate  V.  and  V by 

_ (x)  q^ 


V.  = V + -T 

I HTte  a 


(6) 


Also,  the  interstitial  potential  difference  AV.  is 


AV. 


8a[q.|.-(!rx)qg]+23  (x)qg 


kite  a 
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where  0=0.28/8  for  the  zincblende  structure.^®  In  the  above  expressions 


6 = li(6,+62) 


B = JiCBj-Bj) 


{8a) 

(8b) 
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The  functions  3^{x)  and  6 (x)  are  graphed  in  Figure  2. 

We  may  now  apply  this  model  to  the  calculation  of  the  interface  dipole 

at  a heterojunction.  We  consider  an  ideally  abrupt  heterojunction  between 

compounds  AB  and  CD.  The  parameters  of  the  charge  distribution  depend  only 

on  the  atomic  species  connected  by  each  bond.  For  example,  q (AB),  q_(AB), 

D T 

and  x(AB)  are  the  parameters  of  the  A-B  bonds  in  the  bulk  crystal,  while  qg(CB), 
q.|.(CB),  and  x(CB)  describe  the  C-B  bonds  occur ing  at  the  interface.  To 
evaluate  the  dipole,  we  average  the  charge  distribution  over  the  planes  parallel 
to  the  interface.  The  dipole  calculation  then  becomes  a simple  one  dimensional 
electrostatic  problem.  Explicit  calculations  for  heterojunctions  on  the  (100), 
(111),  and  (110)  planes  give  the  following  results.  The  dipoles,  and  therefore 
the  band  lineups,  are  equal  for  (100)  and  (111)  junctions  of  similar  polarity. 
If,  for  example,  the  B and  C atomic  planes  are  adjacent  to  the  interface,  the 
difference  in  mean  potentials  is  given  by 


AV(lll)=j^  {2q.|.(CB)-q.|.(AB)-q.j.(CD) 
o 

+x(CD)  [l-x(CD)]qg(CD)-x(AB)  [ 1 -x (AB) ] q^ (AB) } 


(9) 


The  first  three  terms  in  the  above  expression  represent  the  transfer  of 
charge  which  actually  gives  rise  to  the  dipole  while  the  last  two  terms 
are  due  to  the  effect  of  the  bond  charge  on  the  mean  potential  in  the  bulk. 
If  the  heterojunction  is  of  the  opposite  polarity,  that  is,  if  the  A and  D 
atoms  are  bonded  across  the  interface,  the  mean  potential  shift  is 
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(10) 


av(TTT)  = ^ 

o 

+x(CD)  [l-x(CD)]qg(CD)  - x(AB)  [l-x(AB) ]qg(AB) . 

For  a heterojunction  on  a nonpolar  (110)  plane,  AV  is  given  by 

AV(llO)  = {q^(CB)-q^(AD)+x(CD)  [l -x(CD)  ]qg(CD) 

o 

(11) 

-x(AB)  [l-x(AB)]qg(AB)} 


If  we  take  the  differences  between  the  above  expressions,  we  find 


AV(1  11)-AV(1  10)  = AV(llO)  - AV(T'i'l)  = 

2^  {q^(AD)  + q^{CB)  - q^(AB)  - q^(CD)} 


(12) 


Therefore,  the  dipole  for  a (llO)  heterojunction  is  equal  to  the  mean  of  the 
dipoles  for  (111)  and  (ill)  heterojunctions. 

An  examination  of  equation  (12)  shows  that  if  q.|.  can  be  written  as 
a difference  between  ionic  parameters,  there  will  be  no  orientation 
dependence  of  the  band  lineups.  This  is  similar  to  the  results  derived 
in  Reference  1.  Also,  for  a heterojunction  between  compounds  which  share 
a common  element,  such  as  GaAs  - AlA*,,  there  will  be  no  orientation 
dependence.  In  such  a heterojunction  on  a polar  orientation,  the  position  of 
the  interface  plane  is  ambiguous,  so  there  is  no  distinction  between  (100) 
and  (100)  heterojunctions,  for  example. 

The  above  expressions  can  be  applied  to  the  prediction  of  heterojunction 
lineups  given  the  parameters  of  the  charge  distribution.  Unfortunately, 
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(13) 


when  we  included  the  bond  charges  in  the  v.ilence  charge  model,  we  also 
included  the  uncertainty  in  the  magnitude  of  those  charges,  and  the  other 
model  parameters  are  rather  sensitive  to  such  uncertainties.  A natural 
way  to  determine  the  model  parameters  would  be  the  following.  First  the 
bond  charge  location  x and  magnitude  q„  should  be  determined  for  the 

D 

bulk  semiconductors  from  the  valence  charge  distribution  calculations. 

The  transferred  charge  q.j.  could  then  be  calculated  from  the  interstitial 

potential  difference  AV.  by  solving  equation  (7).  Chemical  trends  in 

the  parameters  so  determined,  such  as  those  discussed  in  Reference  10,  could 


then  be  used  to  estimate  the  parameters  of  bonds  not  occuring  in  the  bulk 
semiconductors  (such  as  Ge-Ga).  The  sensitivity  of  this  procedure  arises 
from  the  solution  of  equation  (7),  where  small  changes  in  q„  produce  larger 

D 

changes  in  the  calculated  values  of  q.|.. 

In  spite  of  the  difficulties  we  can  make  some  general  statements 
about  the  effect  of  the  bond  charges  on  the  energy  band  lineup.  Only  the 
bond  charges  of  the  bulk  semiconductor  appear  in  equations  (8)  and  (9)"(ll). 

We  may  interpret  those  terms  containing  the  bond  charges  raising  the  bulk 
energy  bands  by  an  amount 

B'^(x)qg  x(I-x)qg 

hirca  * It  a (13) 

o 

If  we  use  Walter  and  Cohen's  estimates  of  qg  (Ge)  = O.lAbe  and  qg 
I (GaAs)  = 0.080^  along  with  the  values  x(Ge)  = 0.50  and  x(GaAs)  = 0.37  from  our 

^ I calculations,  we  get  AEgj.(Ge)  = 0.92  eV  and  AEg^(GaAs)  = 0.1*6  eV.  Therefore, 

? including  the  bond  charge  can  have  a significant  effect  on  the  energy  band 
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lineup.  In  general,  since  the  magnitude  ''f  the  bond  charge  decreases  with 
increasing  ionicity,  we  would  expect  the  bond  charge  effects  to  raise  the 
bands  of  the  group  IV  semiconductors  with  respect  to  the  IIJ~V  compounds  and 
raise  the  lll-V's  with  respect  to  the  l|-VI's. 

III.  Discussion  and  Conclusions 

We  have  presented  a model  of  the  microscopic  charge  distribution  near  a 
semiconductor  heterojunction  which  includes  the  effects  of  bond  charges.  The 
parameters  of  the  model  are  not  yet  sufficiently  well  determined  that  it 
may  be  used  to  predict  actual  heterojunction  band  lineups.  However,  the 
model  does  provide  a basis  for  comparing  the  energy  band  lineups  of 
heterojunctions  of  the  same  chemical  composition  but  fabricated  on  different 
crystallographic  orientations.  We  find  that  heterojunctions  on  polar  (100) 
and  (111)  planes  of  the  same  polarity  should  have  identical  band  lineups. 
Nonpolar  (llO)  heterojunctions  should  have  a band  lineup  equal  to  the  mean  of 
the  (111)  and  (111)  lineups.  Junctions  involving  a common  ion,  such  as  GaAs- 
AlAs,  should  exhibit  no  orientation  dependence. 

A very  interesting  feature  of  the  model  is  the  form  of  the  expression  for 
the  difference  between  interface  dipoles  on  different  orientations  (equation 
(2)).  This  expression  demonstrates  the  role  of  systematic  chemical  behavior 
in  heterojunction  lineups.  For  example,  if  we  could  write 


q.^(AB)  = X(A)  - X(B),  (l4) 


where  X is  a property  of  the  atomic  species  only,  then  we  would  have 


(15) 


Thus,  equation  (12)  is  an  expression  of  t'i<e  deviation  from  chemically 
systematic  behavior  of  the  transferred  charge  q^.  We  may  therefore 
expect  the  orientation  dependence  of  the  band  lineup  to  be  small. 

The  only  experimental  evidence  of  an  orientation  dependence  are  the  old 
Ge-GaAs  data  of  Fang  and  Howard. They  found  a difference  of  about  0.2  eV 
in  the  mean  conduction  band  discontinuity  between  the  different  crystallo" 
graphic  orientations.  However,  the  scatter  in  their  data  for  any  given 
orientation  was  of  about  the  same  magnitude.  It  is  interesting  to  note 
that  if  equation  (15)  were  violated  to  the  extent  of  0.01  electron  (an  estimate 
based  on  the  results  in  Reference  10),  the  heterojunction  lineups  of  the 
different  orientations  would  differ  by  about  0.15  eV.  However,  Fang  and  Howard 
found  the  conduction  band  discontinuity  of  both  the  (ill)  and  (111) 
heterojunctions  to  be  larger  than  that  for  the  (110)  junction,  in  contrast 
to  our  prediction  that  the  (llO)  discontinuity  should  lie  between  the  dis- 
continuities for  the  polar  orientations. 

Finally,  we  should  note  that  q.|.  is  a measure  of  the  net  dipole  moment  of 
the  charge  distribution  associated  with  a band.  Such  a quantity  can  be  defined 
for  any  charge  distribution,  not  just  for  the  model  studied  here.  We  therefore 
believe  that  the  result  given  by  equation  (12),  expressing  the  relation  between 
the  band  lineup  and  the  systematic  behavior  of  the  bond  properties,  would  hold 
for  any  model  in  which  the  bonds  contribute  independently  to  the  charge  distribution. 
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Figure  Captions 


Figure  1.  A schematic  representation  of  the  assumed  form  of  the 
charge  distribution  associated  with  each  chemical  bond. 


Figure  2.  The  Madeiung-type  functions  6 and  6 , which  describe 
the  effect  of  the  bond  charges  on  the  interstitial 
potentials,  as  a function  of  the  bond  charge  position  x. 


